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The interplay between global constraints and local material properties ol chain molecules is a 
subject of emerging interest. Studies of molecules that are intrinsically chiral, such as double- 
stranded DNA, is one example. Their properties generally depend on the local geometry, i.e. on 
curvature and torsion, yet the paths of closed molecules are globally restricted by topology. Molecules 
that fulfill a twist neutrality condition, a zero sum rule for the incremental change in the rate of 
winding along the curve, will behave neutrally to strain. This has implications for plasmids. For 
small circular microDNAs it follows that there must exist a minimum length for these to be double- 
stranded. It also follows that all microDNAs longer than the minimum length must be concave. 
This counterintuitive result is consistent with the kink-like appearance which has been observed for 
circular DNA. A prediction for the total negative curvature of a circular microDNA is given as a 
function of its length. 
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Here, a simple topological sum rule for closed space curves will be applied to a physical model where the molecular 
properties arise by local interactions. The motivation is to better understand what determines the shape of closed 
loops of chain molecules, e.g. DNA. Intuitively, the double-stranded helical geometry of DNA reveals that topological 
restrictions apply and suggests that one can expect the existence of fascinating phenomena. One example hereof was 
recently reported in a study of the supercoiled DNA that forms plectonemes, where long-ranged dynamical effects 
involving two distant plectonemes were observed [T]. Theoretical aspects of the physics of polymers and DNA have 
been reviewed in detail [H[3]. An early triumph is the White's theorem, also associated with Calugareanu and Fuller 
[1HZ] , that imposes powerful restrictions for the ways a closed ribbon can change geometry, and ways it cannot change 
geometry. The theorem describes the relationship between linking, twist and writhe, i.e. Lk = Tw + Wr. The 
White's theorem does not say anything about the material properties of DNA. For example, it does not say whether 
the two strands prefer to be straight, or in a helical configuration, relative to each other. The DNA has a molecular 
preference for the right-handed double helix (B-DNA). This preference breaks the symmetry and makes DNA chiral, 
which gives it special properties. For example, a coupling between strain and twist that we will assume to be local as 
it depends on short-ranged molecular forces. Such properties therefore depend on the local geometry. In most cases 
global topological requirements and local geometrical preferences will be in conflict. In the following we will focus on 
the case where they are concurrent. 

For two three-dimensional space curves, a and /3, that are closed and oriented their Gauss linking number, Lk(a, 0), 
is a topological conserved integer [5] . The linking number describes the number of times the two curves interlink each 
other. The self-linking number, Lk(a,N), of a regular space curve, a, requires a differentiable field of unit normal 
vectors, N. In the following we assume the curve to be smooth to the order of at least C 3 . Together, the unit tangent, 
the unit normal, and the unit binormal vectors frames the curve. The White's theorem applies equally well to framed 
curves and the total twist becomes the twist integral of the vector field N besides a factor of 2tt. Another number 
which describes the vector field N is the winding number. The winding number requires to be measured relative to 
something, i.e., we must have two vector fields defined at the same time. Here we will use the Frenet frame Np, 
actually a modified Frenet frame, N m p, [9 ], as reference. It later becomes a physically relevant frame for winding 
comparisons of intrinsically chiral molecules with non-chiral molecules. 

Let 9 be the angle between the two local frames, i.e. 8 = arccos(A • N m p), then the winding number, $, of N 
relative to the modified Frenet frame becomes 



* = ±f 9 ^s. (1) 
2ir J a ds W 

One could have chosen another reference frame, for example the " Dcnnis-Hannay-Maddocks" frame [TO]. With this 
type of reference frame the winding number of the normal vector field becomes equal to its linking number. With the 
modified Frenet frame as reference, the winding number becomes 



$ = Lk(a, N) - Lk(a, N mF ) . (2) 

Therefore $ is a topological conserved integer. Let us now consider a bundle of frames, N(u), u G R, which are 
differentiable in the vicinity of u = and for which N(0) = N. It follows from differentiation of Eq. |lj that 

/^>d s = 0. (3) 
J a duds 

This is a remarkable simple result that we denote the twist neutrality lemma. In words, the lemma states that the 
differential changes in the rate of the winding progression must obey a zero sum rule. 

Double-stranded DNA is a helical molecule and it therefore has a chiral symmetry, hence some of the properties 
of DNA will reflect this symmetry, e.g. its response to strain. For comparisons of the properties of chiral DNA to 
a hypothetical DNA with no intrinsic chirality it will be useful to choose for the space curve the central line of the 
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double-stranded DNA. Then the framing can be described by the helical path of one of the two strands. By placing a 
ribbon containing the two strands and the central lines of the helical DNA one can see that the linking of the central 
line with one of the strands is equal to the linking of one strand with the other [11] . 

In the following we will make the assumption that the strain-twist coupling for a short stretch of DNA only depends 
on its local geometry, i.e. on its strain, curvature and torsion (cr, k,t). Hence, at small strains the incremental change 
in the rate of winding, 8, will depend only on k and r, 



6 ^ = d^s l -°- (4) 

The physics which underlines this assumption is the relatively local nature of the interactions within DNA between 
its constituent atoms. After all it is their interactions that determine the molecular structure as well as its reaction 
to strain. In Eq. Q we have chosen not to make the orientation of the double helix, i.e. enter the equation. This 
means that material properties are considered uniform over about a pitch of the double helix which is ~10 base pairs 
(bp) long. From applying the twist neutrality lemma to strain we now have 



/ 8(k, t) ds = . (5) 

JDNA 

This restriction, to be twist-neutral in the sense described, is a severe constraint that significantly reduces the number 
of possible closed space curves that double-stranded DNA can make. They are not only topological restricted by the 
White's theorem, they are also geometrically restricted because of Eqs. Q and (J5J). Thus, strained non-conforming 
paths would force the DNA to depart from Eq. fl5j). In the worst case, it could break one or both strands of the DNA 
molecule, in milder cases it could lead to a local detachments of the two strands from each other, i.e. local melting 
of the DNA. Denaturation bubbles have been observed in atomic force microscopy (AFM) images of circular DNA 
[12L 113) and their closure dynamics has been studied by dynamics simulations |14j . 



The implications of twist neutrality can easily be addressed for closed planar curves such as plasmids deposited 
on a substrate. Long plasmids can display many crossings such as supercoils and plectonemic structures [15] . For 
shorter plasmids cooperative kinking at distant sites has been considered [16] . The simplest cases are planar plasmids 
with no self-crossings, i.e. the centerline of the DNA is topologically equivalent to a circle. Recall that straight DNA 
winds when stretched. It rotates opposite of unwinding, a behavior that sometimes is denoted overwinding |17l I18j . 
In the language of this paper it means that 8 > 0. A straight section of DNA can therefore not satisfy the twist 
neutrality criterion. Well, that is a hypothetical statement as a straight line it is not a closed curve, neither! Can 
a circle of DNA fulfill the twist neutrality criterion? As DNA bends, we will assume that 8 decreases monotonically 
as a function of the curvature describing this bending and that 8 eventually becomes negative. This behavior is also 
observed for tubular double helices [TH], while the opposite behavior, i.e. for 8 to increase monotonically, would not 
be possible as we know that DNA can form closed loops. For a perfect geometrical circle to fulfill the twist neutrality 
criterion it must have 8 = everywhere, which is only fulfilled for a specific radius Rq. It means that plasmids that 
are truly circular, and not just called circular because they are loops, come in one specific size with the plasmid length 
Lq = 2t:Rq. Consider a plasmid of length L different from Lq. We series expand 8 to include the linear turn, i.e., 

8 = \(\k\ - no) , (6) 

where K is the curvature with sign as conventionally used for two-dimensional curves, kq = I/Rq is the curvature of 
the circle of radius Rq, and A is a constant that determines the coupling between twist and curvature. Since the curve 
is assumed planar, 8 is independent of the torsion, t. The twist neutrality criterion then becomes 

jfV(*)|-/*)~dt = 0, (7) 
where re(t) is the curvature measured along the plasmid, and s(t) the corresponding arc length. Therefore, we obtain 



l L 
\k\cIs = 2ir — 



(8) 



FIG. 1. Examples of non- convex solutions to the twist neutrality condition for (A) L/Lq = 1.50, and (B) L/Lq — 1.46. There 
are no solutions with L/Lo < 1. 



In two dimensions the integral in Eq. |8]) over a closed curve is always greater than or equal to 2ir [3U] . Consequently, 
there are no possible solutions for plasmids with a length shorter than Lq. Shorter plasmids seem to be required to 
be linear, or to melt and become circular single strands, which presumably themselves have yet a smaller lower cutoff 
for their circumference. 

What are the possible shapes of plasmids which satisfy the "first-order" twist neutrality criterion Eq. (JsJ) with a 
length L larger than the smallest plasmid length, L$l The shape must be concave as the integral of the total curvature 
equals 2-7T exactly when the plasmid is convex [20j . To find a possible shape we cannot use a circle or an ellipse as 
they are both convex. One concave curve with a simple mathematical description is the Cassini oval, which is given 
by the polar expression: 



where a is half the distance between the two foci, and b the square root of the product of the distances to the two 
foci. The shape of the Cassini ovals depends on the ratio e = a/b. For e < 1 the shape is a single loop with an oval 
(convex) or peanut (concave) shape, for e = 1 the lemniscate is produced. Figure 1A shows a peanut-shaped Cassini 
curve which is a solution to Eq. Q with L = 1.5Lq. The parts of the Cassini oval where the curvature is negative 
(the progression of the curve is clockwise) are clearly visible. However, the Cassini ovals are not good descriptors 
of circular DNA which generally are more random in appearance. The desire for some amount of randomness was 
demonstrated in an AFM study of the bond correlation function 21j. In this AFM study it was shown that the 
halfway bond correlation was not completely antiparallel, a result which infers some amount of randomness. By a 
simple random procedure choosing 15 points and generating a closed curve with a spline, a somewhat more realistic 
(though not perfect) concave curve is drawn, see Fig. IB. For this curve L = 1.46Lo- The concave nature of closed 
plasmids is also consistent with the existence of the type of kinks observed in circular DNA [251 H3] • 

For the general case of a concave curve which satisfies Eq. (fsh one has, 



The integral is the total negative curvature without sign, where k_ is equal to \k\ when k < 0, otherwise K- is zero. 
This shows that the total negative curvature must be a linear function of the length, L, of circular DNA. A prediction 
that can be experimentally tested by plotting J_ versus L, e.g. in a scatter plot. Likewise, using a corresponding 
definition for k + , one can calculate a total positive curvature, 



r 4 - 2aV cos(2(9) = b 4 



a 
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(9) 




(10) 




(11) 



hence 1+ — I- = 2w and 2+ + J_ = 2-kL/Lq. One way to estimate the value of Lq would be to invoke molecular 
modeling of DNA on a sufficiently large scale to study the influence of curvature. Tubular models offer a way of 
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getting a first attempt of an estimate. Recently, we found that a tubular model of B-DNA has no strain-twist coupling 
for a circular ring structure with a diameter of 2R = 156 A [24]. The corresponding estimate of Lq is about 124 bp. 
Within the tubular model the uncertainty of this result is at most 5-10 %, however, a tubular model may not fully 
capture all numerical aspects of the molecular structure of DNA and the real discrepancy may well be larger. The 
above estimate for the lower bound, Lq, is in reasonable agreement with lengths recently found for small loops of 
double-stranded DNA in a study of extrachromosomal microDNAs [25] ■ They identified tens of thousands of short 
circular DNA in mouse tissue as well mouse and human cell lines with the overall size distribution in the range 80 to 
2000 bp with a strong propensity for lengths peaked at around 150 and 300 bp for the human DNA and at 180, 360, 
540, 720 and 900 bp (plus perhaps two more peaks) for mouse DNA. Circular single-stranded loops were reported 
as being yet shorter. The apparent tendency for an integer multiple of a preferred length can in the perspective of 
the above analysis be suggestive of microDNAs that coil multiple times when in vivo or are forming plectonemic 
conformations. Alternatively, it could also have other origins such as the specific base pair sequence. 

Theoretically, the regime of relatively stiff polymers for which the elastic response prevails has been considered 
[2"rJl |2"7| . At biological temperatures DNA is fairly close to the melting temperature of the double strand; the melting 
temperature depends in detail on such properties as the salt concentration of the aqueous solution and on the sequence 
of the DNA. Therefore, thermal fluctuations in DNA are omnipresent. Fluctuations in DNA mini-plasmids have been 
modeled as a Kirchhoff-Clebsch rod with no strain [28]. Twist neutrality engender differentiable frames for circular 
DNA, which is suggestive of the existence of well-behaved excitations and fluctuations, e.g. breathing-modes. It 
would be interesting to include strain and stretching in detailed studies of the different kinds of modes and their 
stability [5U1 [30] , and to include twist neutrality in twist-stretch elasticity studies of closed loops of chiral filaments 
[31] - It would also be interesting to include higher-order curvature terms in the twist neutrality condition. This 
includes second-order terms which are also found in bending energy and elastica |32j . 



In summary, a twist neutrality perspective has been given on molecules which display a strain-twist coupling 
reflecting the chiral nature of the molecule. For circular plasmids is was shown that all such plasmids except the 
one of minimal length, Lq, must have a concave shape. The predictions of Eq. (10) can be experimentally tested by 
measurements of the total negative curvature. Further, this will allow for a self-consistent estimate of Lq. 
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